home *** CD-ROM | disk | FTP | other *** search
/ Die Ultimative Software-P…i Collection 1996 & 1997 / Die Ultimative Software-Pakete CD-ROM fur Atari Collection 1996 & 1997.iso / a / a_funk / geograph.tos / GEOGRAPH / UTMKOORD / UTMKOORD.BAS < prev    next >
Encoding:
BASIC Source File  |  1994-02-19  |  1.9 KB  |  34 lines
  1. REM Programm zur Umrechnung von UTM-Koordinaten in geographische Koordinaten
  2. REM z.B. zur Ermittlung des Standorts aus der Deutschen Generalkarte
  3. REM             Christoph Kessler DB1UQ @ DB0GE 19.2.1994
  4. REM Ein UTM-Gitternetz ist vor allem in den "Deutschen Generalkarten 1:200000"
  5. REM zu finden. Damit spart man den Kauf einer detaillierteren Topographischen
  6. REM Karte und das damit verbundene mühsame Peilen mittels Lineal auf die nur
  7. REM am Rand angebrachten Koordinaten. Das Gitternetz besteht aus Quadraten
  8. REM von 10 km Seitenlänge, d.h. 5mm entsprechen 1km, damit ist eine
  9. REM Abschätzung auf 100m (=0,5mm) möglich. Auch eine GPS-Navigations-
  10. REM empfänger ist nicht genauer. Zwei Korrekturfaktoren für die
  11. REM Breitengrade der BRD sind empirisch ermittelt:
  12. REM Zunächst eine lineare Abweichung von 19,73 km
  13. REM Die quadratische Abhängigkeit vom Abstand zum Mittelmeridian
  14. REM korrigiert einen Fehler von ca. 5 km auf den Grenzmeridianen
  15. REM Mit diesen beiden Korrekturen ergeben sich Abweichungen
  16. REM bei ca. 200 ausgemessenen Kartenpunkten von < 500m
  17. REM In dieser Basic-Version ist COS (2*PI)=1 , also nicht in Grad!
  18. REM Hier ist PI vordefiniert, sonst einfügen: PI=3.1415926
  19. REM bitte beides für andere Basic-Dialekte nachprüfen
  20. REM
  21. DEFDBL "A-Z": REM  nur für Atari Omikron.Basic - Double Precision einstellen
  22. MM=9: REM  Mittelmeridian für UTM-Zone 32, östlich vom 12.Längengrad MM=15
  23. GRDPROKM=360/(2* PI *6370): REM         Grad pro km Erdumfang
  24. KORR1=19.73: REM                        Korrektur 1
  25. KORR2=9.359E-5: REM                     Korrektur 2
  26. INPUT "UTM-km zum Äquator (z.B. 5481.5 )          ";UB
  27. INPUT "UTM-km bezgl. Mittelmeridian (z.B. 465.6 )  ";UL
  28. UL=UL-500: REM                          Mittelmeridian entspricht "500km"
  29. GB=GRDPROKM*(UB+KORR1-(KORR2*UL*UL))
  30. GL=MM+(UL*GRDPROKM/ COS(2* PI *GB/360))
  31. PRINT "Geographische Breite:                       ";GB;"  Grad"
  32. PRINT "Geographische Länge:                         ";GL;" Grad"
  33.  
  34.